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We reconsider the proposal of excited dark matter (DM) as an explanation for excess 511 keV 
gamma rays from positrons in the galactic center. We quantitatively compute the cross section for 
DM annihilation to nearby excited states, mediated by exchange of a new light gauge boson with 
off-diagonal couplings to the DM states. In models where both excited states must be heavy enough 
to decay into e + e~ and the ground state, the predicted rate of positron production is never large 
enough to agree with observations, unless one makes extreme assumptions about the local circular 
velocity in the Milky Way, or alternatively if there exists a metastable population of DM states which 
can be excited through a mass gap of less than 650 keV, before decaying into electrons and positrons. 

Dedicated to the memory of Lev Kofman 

PACS numbers: 98.80.Cq, 98.70.Rc, 95.35.+d, 12.60Cn 



I. INTRODUCTION 

There is presently much discussion about whether sev- 
eral anomalous observations in galactic gamma ray and 
cosmic ray astronomy are better explained by dark mat- 
ter models or by more conventional astrophysics. One 
such example is the 511 keV gamma ray excess from the 
galactic center, which has been observed for 40 years, 
most recently by the SPI spectrometer on the INTE- 
GRAL satellite [l[ ; for a recent review see ref. Q . Pulsars 
[ij[ , gamma ray bursts Q , supernovae 0] , low- mass x-ray 
binaries [f|, the galactic black hole [3] and anisotropic 
propagation effects [8| have been suggested as sources of 
positrons whose annihilation could explain the observa- 
tions, but there is no consensus within the astrophysical 
community as to which of these might be the right ex- 
planation. 

There have also been many attempts to explain the 
511 keV signal using particle physics models, including 
annihilations of MeV scale dark matter (DM) [9( or mil- 
licharged DM [ioj . emission from cosmic strings [TTI j . 
decays of sterile neutrinos [ l 2| . axinos [l3| . moduli (or 
modulinos) [Hj], WIMPs [l5 l |lrj| lig ht photinos [13, or 
emission from composite objects [lEj. Excited dark mat- 
ter (XDM) [l9| is a particularly appealing example, in 
which heavy DM particles in their ground state scatter 
into excited states, XoXo — > XiXi) with mass difference 
SM = Mi — Mo. The excited states subsequently decay 
into e + e~xo, where the leptons are approximately non- 
relativistic due to the mass difference being close to the 
threshold value AM > 2m e . If the different DM states 
are members of a multiplet of a new hidden gauge sym- 
metry, this framework has the advantage of being able 
to explain the small mass splittings naturally through 
quantum loop effects [20j; moreover the same class of 
theories can potentially explain excess positrons seen by 
higher energy experiments including ATIC [21], PPB- 
BETS [lljLPAMELA 23] and the Fermi Large Area 
Telescope [24j . 



The XDM proposal requires that the excitation cross 
section be large, in fact close to the unitarity bound, in 
at least a few (and possibly many) partial waves [15j . 
[25| . Ref. [Til made a first attempt to achieve such large 
values in a model where the excitation was mediated by 
exchange of a light scalar <f> with in the MeV— GeV 
range. It was recognized that one must resum ladder dia- 
grams with multiple <f> exchanges when the DM particles 
are scattering at low velocity, but a quantitatively reli- 
able way of doing so was not yet appreciated at the time 
of this work. 

Subsequently ref. [20|] showed how the calculation can 
be set up in the framework of nonrelativistic quantum 
mechanics of a two-state system. It provides an example 
of the Sommerfeld enhancement [2(| [27[, which in the 
last few years has been widely studied in the context of 
galactic dark matter annihilations. However, unlike the 
simpler one-state system where the enhancement can be 
approximated analytically, in the two-state system no an- 
alytic results for the excitation cross section are known 
(however see ref. [27j for recent progress in the annihi- 
lation cross section for multi-state DM). Ref. 28] made 
a first attempt to numerically solve the two-state sys- 
tem and to perform a preliminary scan of the parameter 
space. Technical challenges limited that analysis to rel- 
atively small coupling strengths of the exchanged boson 
to the dark matter. One goal in the present work is to 
extend these results to stronger couplings and to explore 
more widely the range of possibilities, including the de- 
pendence on the mass of the exchanged boson. A second 
is to explore the dependence of the predicted rate on pa- 
rameters of the dark matter density and velocity profiles 
in the galaxy. 

In section|TT]we review the nonrelativistic quantum me- 
chanical formulation of the problem and give a classical 
estimate of the number of partial waves that can be ex- 
pected to significantly contribute to the excitation cross 
section. In section IIIII we discuss the numerical method 
for computing partial wave amplitudes /; as a function 
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of velocity, and present sample results for /; , as well as 
a survey of the range of boost factors which arise from a 
broad exploration of parameter space. In section llVl we 
show how these results go into the calculation of the rate 
of positron production in the galactic bulge; our choice 
of different DM density and velocity dispersion profiles 
is explained there. Section IVl presents the results of our 
scan of the XDM model's parameter space, consisting of 
the interaction strength and the mass of the exchanged 
particle, and the mass splitting between ground and ex- 
cited state. We at first hold the DM mass M fixed at 
a TeV [29j], but then explore the dependence on Mo and 
galactic DM distribution parameters at some optimal val- 
ues of the microscopic parameters. 

Our results show that a large enough rate to match 
observations cannot be achieved in the most straightfor- 
ward implementation of the XDM scenario, but in section 
I VII we review a modified version which can overcome this 
deficit. In this version of XDM, dark matter has at least 
three states; the middle state is stable or metastable and 
can excite over a smaller gap to the highest state. In 
section I VIII we compare the predictions for the angular 
distribution of 511 keV gamma rays with the observed 
signal. Conclusions are given in section rVlIII Appendix 
[^derives the Born approximation for the excitation cross 
section, and IB"! reviews the derivation of the radial depen- 
dence of the DM escape velocity. 



II. EXCITATION CROSS SECTION 

A key realization in ref. [2(| is that the DM particles 
in the galaxy are highly nonrelativistic, and so it is not 
necessary to use the full apparatus of quantum field the- 
ory to analyze their scattering. Instead one can use the 
low energy effective theory, which is quantum mechanics. 
This enormously simplifies the problem, getting around 
the need to resum a perturbative expansion in ladder di- 
agrams to find an effect which is nonperturbative, and 
enhanced by 1/v (or 1/v 2 ), the inverse velocity of the 
DM particles. 



A. Quantum mechanical analysis 

We assume that the DM is Majorana and that its mul- 
tiplicity arises from it being in a nontrivial representation 
of a nonabelian gauge symmetry. This immediately im- 
plies that any two members of the multiplet, say Xi an d 
Xi-, can only interact via off-diagonal couplings to the 
gauge boson. Furthermore the gauge symmetry should 
be spontaneously broken so that xi an d X2 can be non- 
degenerate; then the exchanged vector boson gets a mass 
\x. This leads to an attractive Yukawa potential, but only 
as an off-diagonal term in the matrix elements of the po- 
tential for the two-state system. It is important to note 
that each state consists of two particles, 
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|Xi>Xi), 



|X2, Xi) 



by virtue of the fact that whenever a gauge boson is ex- 
changed, the color of both DM particles has to change. 
There also exist the states |xi, X2) and \xi, Xi) which live 
in their own superselection sector, but we do not need to 
consider them because xi is presumed to not be present 
in the initial state. The diagonal terms have no interac- 
tions; this part of the Hamiltonian consists only of the 
the mass-energies of the two particles. Subtracting the 
mass of the ground state, 2Mq, the matrix potential in 
the basis |T]) is 



Vij = 







-a a e 



-fir 



-a g e tir jv 



2SM 



It 



(2) 



where a g — g /4ir is the fine structure constant for the 
SU(N) gauge coupling g. The wave function for the two- 
state system (with components labeled by index i) in 
the CM frame is = Y,i p i ( cos0 ) R ki( r )^ where k is 
the initial momentum. Defining 3>i,j(r) = UjL/r, the 
Schrodinger equation is 
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M ~ M ' V M or 2 "' J ' "7 "" J M 

For the numerical solution it is useful to rescale r = 
(a g /28M)x and define the dimensionless variables 



1(1 + 1) 
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r = Mi 
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where v t — y / 25M /Mo is the threshold velocity for pro- 
ducing the excited state. Then the Schrodinger equation 
takes the form 

/ 1(1 + 1) 



V = 



r(v - A) $ = 



1 



(5) 



(6) 



with the dimensionless potential V. We require A > 1 
for the initial state to have enough energy to produce the 
heavier |x2,X2) final state. 

To extract the scattering amplitudes, we decompose 
the numerical solution into incoming and outgoing waves, 
<&°" t and Partial wave unitarity implies the 

conservation of flux, /c|$^| 2 = k\$ff\ 2 +k f \$?$ > \ 2 (where 
k' 2 = k 2 — 2Mo6M), which we use as a check on our 
numerics. The fraction of incoming |xi,Xi) states which 
gets converted to the \x2, X2) Anal state is thus 
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(7) 



in the Ith partial wave. // provides a measure of the 
extent to which a given partial wave can saturate the 
unitarity bound // < 1 for 07, its contribution to the 
cross section. The partial wave cross section is 



n(2l + l) 
-1 = -Mf^ fliv) 



(8) 



(1) and the total cross section is a — <7j. 
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B. Classical treatment 

To get additional insight, it is useful to think about 
the classical version of the problem, in the limit where 
the scattering is elastic. The standard method for solv- 
ing the central potential problem is to change variables 
to it = 1/r and to solve for 9(u) instead of r(t), where 9 
is the polar an gle in the scattering plane (see for example 
section 4.5 of [30|). Take M — \Mq to be the reduced 
mass. Then it is straightforward to solve the first inte- 
gral of the motion for 9(u) with initial conditions where 
the particles have velocity v = v Te \/2 in the center-of- 
mass frame and impact parameter b. Taking w = bu, the 
scattering angle is 



dw 1 



2a 



Mbvi 



3 e ~fJ.b/w 



-1/2 



(9) 



where wq ~ 1 is the turning point of the effective po- 
tential, the point where the integrand diverges. Notice 
that 9 S = in the limit where the scattering potential 
vanishes. The result © shows the classical origin of the 
Sommcrfcld enhancement, where the effect of the poten- 
tial is strengthened by the factor 1/v 2 at low velocities, 
since the particle has more time to be influenced when it 
is moving slowly. 

We can estimate the maximum angular momentum 
I = Mgvb which gives significant scattering by demanding 
that (2a g / Mbv^ cl )e~ fJ ' b >, 1. Solving for the argument of 
the exponential gives 
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(max ~ V In 
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The prefactor M^v/ji is what one would have obtained 
by assuming the range of the force is Substituting 
this estimate for / max (without the logarithm) into the 
formula for the cross section © and assuming fi = 1 for 
all I < l max results in the geometrical value of the cross 
section, a = n//j 2 . The extra In (a ff /x/Mof 2 ) factor is 
the origin of the Sommerfeld enhancement, which boosts 
the cross section by the square of the logarithm for low- 
velocity scattering. 

For the inelastic process, these classical insights need 
to be modified, since there is a threshold v > v t for pro- 
duction of the excited state. Nevertheless they provide 
some idea as to how many partial waves one may expect 
to be important in the quantum mechanical cross section 



III. PARTIAL WAVE AMPLITUDES 

In this section we describe our method of numerical so- 
lution of the Schrodinger equation, leading to the partial 
wave amplitudes fi(v). In each partial wave, one bound- 
ary condition is that there are only incoming waves of 
state |1) from r — oo, whereas the outgoing waves are 



an admixture of |1) and |2). Specifically, |2) must ap- 
proach an eigenstate of radial momentum with positive 
eigenvalue, 



O r 



1.2 



(11) 



where k' = ^/k 2 — 2Mq5M . This constitutes two con- 
ditions, one for the real part and one for the imaginary 
part of the wave function. At the origin, both compo- 
nents behave like r l+1 , but their relative amplitudes are 
not known, so one must parametrize ~ rl+1 {\) with 
b a complex number. We thus have a shooting problem: 
the real and imaginary parts of b must be adjusted so as 
to satisfy the two boundary conditions at infinity. Stan- 
dard algorithms exist for this kind of problem; we adopt 
the routines in ref. 
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Instead of shooting, one can alternatively use a sim- 
pler approach, which is to solve the Schrodinger equation 
with the correct r l+1 behavior near r — 0, but arbitrary 
relative amplitudes of and $1,2, to obtain some solu- 
tion $o, which does not have the right behavior at large 
r. Since the complex conjugate $q is also a solution, 
one can algebraically construct the b linear combination 
a$ + b$Q that has the desired behavior at large r. This 
is much faster than shooting because no iteration is re- 
quired. 

In principle, the method is straightforward, but com- 
plications arise when one tries to consider parameters in 
the regime Y > 1, that is, a g > y/2SM/M . Notice 
that this need not be a particularly strong coupling for 
\ the application we have in mind, where SM ~ 2m e and 

J (1°) M - TeV: then a g > 1CT 3 . The algorithms break down 
except when 77, which determines the dark gauge boson 
mass, is sufficiently large. The problem arises from the 
need to consider very different scales 1/r and I/77; in 
particular the wave function oscillates many times in the 
interaction region x < 1/77. Then (we suspect) the solu- 
tion with the desired behavior at large r becomes an ex- 
ponentially small component of the generic solution, and 
so its extraction gets lost in the numerical noise. Despite 
this limitation, we will be able to explore the parameter 
space widely enough to see how the results extrapolate 
to the difficult regions, and thus give a complete charac- 
terization of the solutions. 



A. Results for /; 

We show two examples of the results for eq. ([7J in 
figure [U as a function of A — 1 = v 2 /v 2 — 1; recall that 
v is the velocity of one of the particles in the center of 
mass frame, and v t is the threshold velocity for producing 
the excited DM states. These results look quite different 
from the usual single-state s-wave Sommerfeld enhance- 
ment, which has numerous sharp resonance bands at low 
velocity. Here, there can be at best a few oscillations 
in a given partial wave, but not strongly peaked. These 
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FIG. 1: Top: fraction of partial waves converted to the excited state as a function of center-of-mass velocity (relative to the 
threshold velocity vt) for two sample sets of parameters, F = 10, r) = 1.2, and T — 340, r\ = 8.2. Each curve is labeled by its 
partial wave Bottom: same, but showing region of smaller velocities. 



might be due to having an approximately integral num- 
ber of wavelengths within the width of the potential well, 
which of course is not square and so one would not ex- 
pect nearly bound resonances to be very sharp. It would 
be interesting to find an analytic approximation to help 
better understand these results. We tried to develop the 
WKB approximation in this context, but without suc- 
cess. 1 

For clarity in figure Q] we have chosen examples where 
the number of partial waves which contribute signifi- 
cantly are relatively small, but for our actual computa- 
tions, there are examples (particularly for small values of 
rj <C VT) where hundreds of partial waves are relevant. 
Using the estimate (flQ|). and the definitions (j4|), we find 



(12) 



We typically computed up to l max = 500 to insure con- 
vergence of the sum over / in the cross section, and we 
restrict ourselves to gauge boson masses, hence values 




of r\ ~ VT, which are large enough to justify neglecting 
higher values of / (it also turns out that the interest- 
ing functional dependence of a on -q occurs in the region 
where rj ~ %/T). It will be apparent from the results 
that this is a mild restriction, in terms of understanding 
trends across the full parameter space of the model. 



B. Boost factors 

There is a large literature on the Sommerfeld enhance- 
ment for s-wave annihilation of DM at low velocities, 
where the effect is characterized by a boost factor, defined 
as the ratio of the enhanced cross section to the Born ap- 
proximation value. Even though all our results for the 
scattering cross section are computed directly, without 
reference to the perturbative value, one might neverthe- 
less be interested to know their ratio in the present con- 
text. To make the comparison, we need the perturbative 
expression for a (here computed using relativistic quan- 
tum field theory), 



1 At least in some regions of parameter space, the WKB approx- 
imation fails due to the breakdown of matching conditions at 
the classical turning point. It is possible that more sophisti- 
cated approximate analytical techniques such as those used in 
[27ll ameliorate this problem. 
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2A 1 (M 2 « t 2 (2A- 1) + ji 2 )v 
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2^ 2 M 2 v 2 (2A- 1) 



(13) 
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FIG. 2: Boost factors as a function of the dimensionless squared velocity A = v 2 /vf — 1 (recall that vt is the threshold velocity) 
for a wide range of the parameters denned in eq. @, T — 10, . . . , 10 5 (rows) and 77 ~ few Xs/T (columns). We assume Mo = 1 
TeV, and the curves in each panel correspond to 8M = 0.2, . . . , 1 MeV, from top to bottom. 



where e = 2M 2 ti^A(A - 1)/( M 2 + M 2 v t 2 (2A - 1)), and 
we have assumed SM <C \x <S Mo to simplify the expres- 
sion; see appendix [A~l for details. 

In figure [2] we plot the resulting boost factors as a 
function of center-of-mass velocity of the incoming par- 
ticles, for the same range of parameters that we explore 
in the next section as being relevant for positron pro- 
duction, namely T = 10, . . . , 10 5 and r\ ~ fewx\/T, for 
SM = 0.2, 0.4, . . . , 1 MeV, and assuming a DM mass 
of Mq = 1 TeV. From the figure it is apparent that the 
"boost factor" is an enhancement only for the smaller val- 
ues of r, and is actually a suppression factor for T > 100, 
in the range of velocities where we have computed. One 
can see the onset of resonant enhancement at low ve- 
locity in the cases of T = 10 and 100, where the boost 
factor suddenly increases as a function of the dimension- 
less gauge boson mass parameter rj. For fixed parameters 
r, rj, A, we see from each family of curves that smallest 
values of SM give the largest cross section. This can be 
understood as being due to Sommerfeld enhancement, 
since for fixed A decreasing SM corresponds to decreas- 
ing the particle velocities. 



One check on these results is provided by considering 
when the Born approximation for our potential should 
be valid. By demanding that the scattered wave be small 
compared to the incoming wave in nonrelativistic quan- 
tum mechanics, and working in the limit of low velocity, 
one obtains the constraint 2Mq\ J q °° rV(r)dr\ <C 1, which 
implies 2a g Mo <C /it, or in terms of the parameters ((4]), 
2r -C 77. This is violated everywhere in the region of 
parameter space we have considered, so there is no con- 
tradiction between our results and expectations based on 
the Born approximation. 

IV. RATE OF POSITRON PRODUCTION 

The rate of positron production within a radius r c of 
the galactic center is given by 

1 f r " 

^e+=x/ (<7v Icl )n 2 (r)4irr 2 dr (14) 
1 Jo 

where n(r) — p(r)/M is the DM number density and 
p(r) is its mass density. The factor of \ is to avoid 
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double-counting, and we have assumed that the DM is 
distributed with spherical symmetry. Since the INTE- 
GRAL signal has a full-width at half-maximum of 8° 
we take r c to be the radius which subtends 8° at 
our distance of Rq — 8 kpc from the galactic center, giv- 
ing r c — 1.1 kpc. We will discuss the uncertainty in Rq 
and its impact on our results below. The predicted rate 
(P3J is to be compared to the observed one (see ref. [l9| 
for a detailed discussion) of i? bs = 3.4 x 10 42 for the 
component coming from the galactic bulge [l9j . 

The cross section in (fT4"f is averaged over the veloci- 
ties of the DM particles, using the Maxwell-Boltzmann 
distribution function 



f(v,r)=N 



v 2 exp (— v 2 /2a v (r) 2 ) , v < v c 



0, 



V > Vr. 



(15) 



with a cutoff at the escape velocity for galactic DM. The 
normalization factor is given by 



N-\t) = 2 5/2 7ra.i!(r) f^erf(A) - Xc 



-A" 



(16) 



where X = v csc /(V2a v (r)). We follow ref. |T3| m us- 
ing the r-dependent escape velocity (see appendix [B] for 
derivation) , 



1 - In 



r < r_ 2 
r > r- 2 - 



(17) 



v c is the circular velocity of DM, assumed to be approxi- 
mately constant at radii near the characteristic scale r_2, 
which is defined below. We take the fiducial value for v c , 



220km/s 



(18) 



as in ref. (l9( and many other references, but below 
we will also explore the sensitivity of our predictions to 
changes in this value. To average over the cross section, 
one integrates over both DM particle velocities, 



(<rv le i)(r) = 87r 2 / dvi / dv 2 
Jo Jo 

f(vi,r)f(v 2 ,r) av Tcl 



d(cos(9) 



(19) 



where u re i — (vf + v 2 — 2v\v 2 cos 9) 1 / 2 , and the cross sec- 
tion is taken to vanish if ^v vc \ is less than the threshold 
velocity defined in eq. ((4]). We perform the integrals nu- 
merically. 



Einasto, a = 0.17 

cored BDM 

■-■ NFW 




0.4 0.6 

log 1Q r/kpc 

FIG. 3: NFW, Einasto, and cored BDM density profiles, in 
the inner 1—10 kpc. Inset shows the region 0.01 — 1 kpc. 



(PDM), and others which contain a baryonic component 
in addition to the DM (BDM). The PDM Aquarius sim- 
ulation [33| finds the relation a 3 oc r 1-875 p(r), where p is 
the DM density. Thus one needs to specify the density 
profile in order to fix a v (r). 

We consider two widely-used density profiles, the NFW 
form 



p(r) 



4p-2 

r-2 



r-2 



and the Einasto form 

p{r) = p-2 exp 



- - 

a \\r-2 



1 



(20) 



(21) 



Here r_ 2 is the radius at which the logarithmic slope of 
the density is —2, the normalization is defined by p_2 = 
/o(r_a), and a is the shape parameter. A fit of the NFW 
form to the Milky Way galaxy in ref. [13] gives 2 r_2 = 
14.1 (h/0.7)- 2 / 3 kpc and p_ 2 = 0.13 (h/0.7) 2 GeV/cm 3 
where h is the Hubble parameter. These are close to the 
best-fit values from the Aquarius simulation for the Aq- 
A-l galaxy, which was the highest resolution simulation 
considered, given in tables 1-2 of ref. 33] , 



r-2 = 15 kpc 
a = 0.17 



(Aquarius- A-l parameters) (22) 



A. Pure DM radial distributions 

The r dependence of the velocity dispersion a v (r) and 
also the escape velocity v csc (r) appearing in f| 15[) de- 
pend on the shape of the DM density profile p(r). We 
will consider several hypotheses for the form of p(r), in- 
ferred from A-body simulations, some based on pure DM 



2 In the notation of ref. I'{4| . r_2 = r s = r v /c and 4p_2 = & c pc, 
where r v is the virial velocity, c is the concentration parameter, 
S c = 100e 3 /[3(m(l + c) — c/(l + c)] and p c is the present critical 
density of the universe. r v is related to the virial mass M v by 

= / 3A/„ \ 1/3 
Tv ~ \4n ■ lOOpc) 

Using their best fit values c = 18, M v = 9.4 X lO n M (see 
erratum of |34H ) leads to the values quoted for r B and p_2- 
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PDM, Einasto \ 
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FIG. 4: Velocity dispersion profile for dark matter + baryons 
(BDM), solid curve, and for pure dark matter with Einasto 
or NFW density profiles, dashed curves. 



We adopt these as our fiducial values. The resulting ve- 
locity dispersion is then 



r-2 



0.625 / / sx 1/3 



P-2 



(23) 



where the velocity scale vo (called cr max in Table 1 of ref. 
331) is determined to be 



v = 260km/s 



(24) 



The radial dependence of a v is plotted in figure HI show- 
ing that the Einasto profile leads to much higher veloc- 
ities in the inner regions r < r_2 than does the NFW 
profile. 



B. Effect of baryons on distributions 

The Aquarius simulation ignored the effects of baryons 
on the evolution of galaxies, but recent studies have in- 
cluded them and shown that their effect is to increase 
the velocity dispersion in the inner region r <, 10 kpc 
(2l-jH3- Fitting the results of figure 2 in ref. we 
find the velocity dispersion 



a v (r) = v 



r-2 



-1/4 



(25) 



where vq happens to take the same value as in Its 
shape is plotted in figure 0] The rise in v toward the 
galactic center can greatly boost the excitation rate of 
DM relative to the pure dark matter case, (|2"3"1) . and we 
will employ both for the purposes of comparison. 

There is some disagreement between the different BDM 
simulations as to the impact of the baryons on the DM 
density profile. Ref. [35j finds that the cusp of the den- 
sity distribution is softened (we refer to it as a "cored" 
profile) , while ref. [37} obtains profiles that are consistent 



with the Einasto form. We will consider both possibili- 
ties in our computations of the positron production rate. 
To study the possible softening effect, we have digitized 
the BDM density profile at redshift z — in fig. 1 of ref. 
[IH and fit its logarithm to a quartic polynomial, 



logio 



71 = 



log^r/kpc) 



(26) 



with coefficients a n = 1.7487, 0.0395, -2.537, 1.459, 
—0.3448, respectively, for n — 0,...,4. This profile is 
plotted along with the NFW and Einasto forms (|2Tj ]) -(|2"T j) 
in figure[3]in the region 1 kpc < r < 10 kpc. For larger r, 
the three shapes are in closer agreement, while for r < 1 
kpc the BDM profile is taken to be constant. 

Other BDM simulations do not find the erasure of the 
cusp; ref. [37J finds good fits to the Einasto profile with a 
range of a values that are consistent with (|2"2"1) . but that 
go as low as a = 0.145. To illustrate the effect of the a 
parameter on the rate, we will consider this lower value 
in addition to the higher one a — 0.17 in (|2"2")l . 



C. Consistency of velocity distributions 



Strictly speaking, the velocity distributions (|23I25[) 
may not be self-consistent when varying the dark mat- 
ter density profile. A better approximation to the self- 
consistent distribution could be derived from the Jeans 
equation, following the technique of [38], for example. 
However, our results are much more sensitive to the shape 
of the density profile than the detailed velocity profile, so 
an improved treatment of the velocity profile should only 
lead to qualitatively small changes in our results. Figure 
[5J for example, shows that the rate of positron production 
is relatively insensitive to the scale of the BDM velocity 
distribution vq near or above the standard value (|24[) . 

We do however, find an enhancement in positron pro- 
duction in going from the PDM velocity dispersion profile 
(f2"3"f to the more cuspy BDM profile (J25J. We therefore 
might worry about overestimating the positron produc- 
tion rate if a self-consistent treatment removes the cusp 
in the profile. We believe it is more conservative to con- 
sider the cuspy profile (l2"5j) . as we are interested in de- 
termining whether or not the XDM mechanism is viable 
for 5M > 2m e . Even with the advantageous profile (|25l) . 
we still find that the positron production rate is only 
sufficient for smaller SM. 



V. SURVEY OF PARAMETER SPACE 

There are four important dimensionless parameters 
that determine the rate of scatterings to produce the 
XDM states in the galaxy. Two have already been de- 
fined in eq. ((U), T and 77, which depend only upon the 
microphysics, i.e., the DM mass, mass splitting, gauge 
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Contours of log 10 R / R obs 
DM + baryons velocity dispersion 

Einasto profile, a = 0. 145 




Contours of log 10 R / R obs 
pure DM velocity dispersion 
. [r 1.875 p(r) ]l/3 
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FIG. 5: Contours of logarithm of predicted over observed rate of positron production in the plane of r/SM for several DM 
velocity and density profiles. (77 is proportional to the exchanged gauge boson mass; see eq. Q.) Those labeled "0" (thickest 
contours) match the observed rate. The DM mass is assumed to be Mo = 1 TeV. 



boson mass, and interaction strength. The others de- motivating us to define 
pend on the ratio between the threshold velocity vt or 

the escape velocity v csc and the parameter vq (that con- 2SM v 2 v 2 

trols the average speed of DM particles in the galaxy), a = — — ^ = b = (27) 

M v£ i>5 
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FIG. 6: Log of predicted over observed rate contours in the plane of mass splitting 5M versus Mo (DM mass), vq (velocity 
dispersion at r_2 = 15 kpc) or v c (circular velocity at r_2). The contour labeled "0" corresponds to the observed rate of 
positron production. The distance to the galactic center is fixed at Ro = 8 kpc in the top right panel (and the left ones), 
whereas it varies linearly with v c in the bottom right one. Other parameters are fixed at F = 10 4 , 77 = 50, a — 0.145 to 
maximize the rate, and the BDM velocity profile is used. 



Dependence on a, b arises from averaging the cross sec- 
tion over the velocity distribution of the DM. In ref. [ll|, 
a simpler estimate of the rate was made using the isother- 
mal distribution function, f(v) — Nv 2 e~ v < v °, giving 

> = #E( 2 '+d f^K^m (28) 

M§vq J v 2 

where v csc is the escape velocity, also taken to be constant 
and in the range ~ 500 — 600 km/s. 3 In the limit of large 
a, b, the rate is exponentially suppressed, ~ (eT a — e~ b ). 
Although we are calculating the rate more quantitatively 
in this paper, the simpler approach makes it clear how to 
expect the results to depend on a, &, at least semiquan- 
titatively. In terms of dimensionful constants and the 
dimensionless ratios, one can parametrize the rate in the 
form 

Re+=^r :1 9(r,V,a,b) (29) 

where the function g of the dimensionless variables has 
complicated dependence on T and 77, but roughly (e~ a — 
e~ b ) dependence on a, b. The contours of log(e~° — e~ b ) 



3 This estimate assumes that v csc 2> vq so that N ~ v 3 . In 
the unphysical limit v eB c *C vo, the normalization factor goes to 



in the a-b plane are roughly linear over some ranges, sug- 
gesting that there should be a quasilinear degeneracy be- 
tween the parameters SM and v^ sc , which we will observe 
below. 



A. Dependence on T, r\, SM 

For our initial exploration of parameter space, we com- 
puted the scattering rate for values T — a 2 Mo/25M = 
10 1 , 10 2 , . . . , 10 5 ; this is the parameter that is most di- 
rectly associated with the strength of the dark gauge in- 
teraction since it multiplies the potential in ©. For each 
value of r, we explore a range of 77 = a g /i/25M (recall 
that 77 determines the range of the interaction since \x 
is the gauge boson mass) that includes the region where 
the rate of excitations is maximized. For each value of 
r, we calculated contours of the positron production rate 
R e + in the 77-a plane. However to make the presentation 
more concrete, instead of using the dimensionless vari- 
able a, we momentarily assume that the model should 
also account for the positron excesses seen by PAMELA 
and Fcrmi/LAT, the latter of which suggests that Mq = 1 
TeV [1^]. This enables us to plot the contours in the 77- 
8M plane. An important question is whether the rate 
can ever be large enough if 5M > 2m e since in the most 
natural models, each excited state must have enough en- 
ergy to decay into an e + -e~ pair. 

In hgure [5] we display results for three cases with the 
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BDM velocity dispersion (|25|) and different assumptions 
for the density profile, and one example with the PDM 
dispersion (f23|) . The largest rates are obtained in the 
former case, using cuspy density profiles, and T = 10 4 . 
Even in this most favorable case however, the mass split- 
ting cannot exceed 600 keV if the predicted rate is to 
match the observed one. The cored BDM density profile 
and the PDM velocity profile give much smaller rates, 
with correspondingly smaller upper limits on 5M, 170 
keV and 200 keV respectively. 



^ 5M-100 keV 
AM > 1 MeV 

FIG. 7: Spectrum of states for inverted mass hierarchy. 



B. Dependence on other parameters 

Ideally one would like to find examples where the pre- 
dicted rate can match observations for mass splittings 
5M >, 2m e , since in the simplest models each excited DM 
state must decay into the ground state plus e + e~ . To- 
ward this end, we have explored the dependence on other 
parameters which we previously held fixed. For this pur- 
pose we fix instead the microscopic parameters at one of 
the highest-rate examples found in fig. [5] Namely we 
take r = 10 4 , r\ = 50, using the BDM velocity pro- 
file, and Einasto density profile with a = 0.145. We 
vary the DM mass M , the velocity dispersion parameter 
yo, and the circular velocity parameter v c - Recently ref. 
[391 ] presented evidence, based on trigonometric paral- 
laxes and proper motions of masers in star-forming re- 
gions of the Milky Way, favoring a surprisingly large 
value v c = 254 ±16 km/s. The claim has been criti- 
cized in subsequent references [40l l4lj ; in particular ref. 
[4~fl | notes that the method used only constrains the ratio 
v c /R = 30 km/s/kpc, where i? is the distance to the 
galactic center. We therefore also consider the variation 
of v c and Ro together. 

The results are shown in figure [5] The top left figure 
shows that there is an optimal value of Mq = 600 GeV 
for maximizing the rate, where SM = 650 keV gives the 
observed rate. This can be understood from the analytic 
approximation (|29p where for fixed SM, the rate scales 
as a 4 e~ a ~ M Q ~ 4 e~ c / M ° . The bottom left figure shows 
that the dependence on vq, which controls the overall size 
of the DM velocity dispersion, is quite weak for values 
above the standard one vq = 220 km/s (see footnote [3] 
for explanation). As expected, decreasing v below this 
value only reduces the rate. 

The right panels of figure[6]show the dependence on the 
circular velocity v c , which controls the escape velocity of 
the DM. Keeping Rq fixed, the upper right panel demon- 
strates that v c would have to be increased to 350 km/s 
to allow for SM as large as 2m e . However increasing Rq 
in a correlated way, as suggested by ref. [4l| , allows for a 
stronger effect, because the cutoff r c on the radial inte- 
gration in eq. (|14[) increases proportionally with Rq. In 
this case, one would need v c = 280 km/s and Ro = 9.4 
kpc. Such a large value of v c is 4cr away from the recent 
mean determination of 236 ±11 km/s [4(| and that of 
Rq is 2.5tr and 3.1er outside the preferred values of refs. 



[H HI respectively. 

We conclude that, by taking all parameters and dis- 
tribution functions to their limits, it may be marginally 
possible to excite dark matter with a mass splitting as 
large as 2m e , although most practitioners would prob- 
ably regard the required values of v c — 280 km/s and 
Ro = 9.4 kpc as being unreasonably large. 



VI. ALTERNATIVE XDM SCENARIOS AND 
TEV DARK MATTER 

In ref. [28| it was pointed out that a long-lived species 
Xi with a mass splitting AM ~ 2m e above the ground 
state xo could be excited to a third state xi with split- 
ting SM <C 2m e in an "inverted mass hierarchy" scenario, 
whose spectrum is shown in fig. [7] Then the XDM mech- 
anism can work via the excitations XiXi ~~ ^ X2X2 followed 
by the decays xi ~^ Xo<s + e~. The reduced value of SM 
allows the rate to be large enough to match observations 
without requiring extreme values for parameters of the 
galactic velocity profile. 

Simple models of nonabelian DM with hidden SU(2) 
gauge symmetry were constructed in ref. (44| to realize 
this possibility. The scenario comes with the expense 
of needing a late-time nonthermal origin for the DM, as 
explained in [28| (see also (45|), since otherwise the same 
process needed for positron production in the galaxy will 
depopulate the xi states in the early universe. If the 
DM goes out of kinetic equilibrium early enough, a small 
remnant of its initial relic density can be maintained; ref. 
(25| optimistically estimates that ~ 1/10 of the initially 
produced relic density can survive. In this case we would 
need parameters corresponding to a rate R e + /R bs = 100 
to compensate for the p 2 ~ 0.01 suppression in the rate, 
relative to the assumed case of a standard relic density. 
The upper right-hand panel of figure [5] (corresponding to 
BDM with a steep Einasto profile) reveals such examples 
(the contours labeled "2") when V ~ 10 3 -10 4 and SM ~ 
50 keV. 

In the case of nonthermally-produced intermediate 
states, in order for this component of the DM to dominate 
over the conventionally produced thermal component of 
the ground state, it is necessary to have an annihilation 
cross section that is stronger than the one needed for the 
standard thermal relic density, which was computed for 
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FIG. 8: Relation between r, f and <j> for the line of sight 
integral p3|) . 

SU(2) DM in ref. This puts a constraint on the dark 
gauge coupling of triplet DM, 

a * > °' 03 (rSv) (30) 

If we assume that M = 1 TeV, and 200 MeV < fi < 
1 GeV, the preferred value for explaining excess leptons 
seen by the PAMELA and Fermi/LAT detectors, this 
puts constraints on the dimensionless parameters T, r/: 

« < r * (ioffffiv) < 5100 < 31 > 
a03 < " * (i5nbv) < 5 < 32 > 

where we have also imposed that a g < 1. Referring to 
fig- El w e see that these can be satisfied for the top panels, 
which use Einasto profile and BDM velocity dispersion, 
even at Y = 10. There thus seems to be considerable 
room in the parameter space for the XDM explanation 
of low-energy galactic center positrons, if one accepts the 
inverted mass hierarchy hypothesis and a relic density of 
intermediate mass states. 



where f is distance along the line of sight, and r 2 (f, (f) — 
Rq + f 2 — 2rR cos (j), as illustrated in fig. [5] The observed 
profile is broadened by the SPI instrumental resolution 
of 3° full width at half maximum, so we smooth it using 
a gaussian, 

L obs ((f>) = N(c) J dfo e-<*-*°?) L(fo) (34) 

where c = 7r(3°/180 o )/(2 3 / 2 In 2). 

Our results for the angular profile are shown in fig. |9l 
for sample microphysics parameters Y = 100 and r/ = 3.5. 
We find that the shapes are relatively insensitive to these 
values, and rather depend more strongly on the DM den- 
sity and velocity profiles, and the parameters a, b defined 
in ([27]) . For the first five panels of this figure, we have ig- 
nored the overall rate and normalized all curves to match 
the observed profile at the galactic center, to better vi- 
sualize the differences in widths of the predicted signals. 
The bottom right panel shows the real rate correspond- 
ing to each assumed density profile, as a function of 8M 
(assuming Mo — 1 TeV as before). 

One can see that for the larger values of 8M, the pure 
dark matter velocity profiles (as well as cored BDM) pro- 
duce too wide a profile to match the observations, while 
the cuspy BDM case gives a narrower result that can 
be consistent if the positrons travel a distance of order 
1 kpc before annihilating, as has been argued is possi- 
ble (46|. This indicates that the kinematic advantage of 
having higher central velocities in the BDM scenario is 
more important than the boost to the cross section that 
one would expect from the Sommerfeld enhancement at 
lower velocities. Thus not only are the cuspy BDM pro- 
files more easily able to match the observed rate, but 
they are also more consistent with the observed angular 
distribution. 



VIII. CONCLUSION 



VII. ANGULAR DISTRIBUTIONS 

Further constraints can be obtained by computing the 
expected angular profiles of the produced positrons. This 
gives one-sided bounds, because there is considerable un- 
certainty in the distance that a positron could propagate 
away from the galactic center before annihilating. But a 
model that predicts too wide an angular profile can be 
excluded since propagation effects will never make the 
true profile more narrow. 

Neglecting propagation effects, the angular distribu- 
tion of positrons is given by an integral similar to that in 
(H3J), but instead of integrating over the volume, one must 
integrate along the line of sight. We define a luminosity 
L(0)by 

L{4>)= / (<jv lci ){r)n 2 (r)dr (33) 
Jo 



In this paper we have tried to give the most quantita- 
tive treatment to date of the excited dark matter mech- 
anism for producing positrons at the galactic center. A 
main technical improvement was the numerical compu- 
tation of the excitation scattering cross section, which is 
nonperturbative and guaranteed to satisfy unitarity con- 
straints. A second difference relative to previous treat- 
ments is that we included the effects of baryons on the 
velocity dispersion of DM in the inner kpc of the galaxy, 
which can significantly boost the rate of excitations. In 
addition we considered NFW and Einasto profiles for the 
DM density, varying parameters determining the cuspi- 
ness of the distributions, as well as those affecting the 
velocity distributions. 

Even making all the most optimistic assumptions for 
increasing the rate of DM excitations, we were not able to 
find realistic parameter values which yield a large enough 
rate, if the mass gap between the ground state \o and 
excited state xi of the DM is sufficient for the subsequent 
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FIG. 9: Comparison of the predicted and observed angular profiles for V = 100, r\ = 3.5 and five values of SM, also varying 
the density profiles. All curves are normalized to 1 at r = to emphasize the differences in shape. Bottom right graph shows 
dependence of the actual rate on 5M for these same density profiles, assuming that Mo = 1 TeV. 



decay xi Xo e + e~ to produce the positron. The largest 
mass gap we found consistent with the observed rate was 
SM = 650 keV, assuming the standard values v c = 220 
km/s and Rq = 8 kpc, respectively, for the circular ve- 
locity of the sun around and its distance to the galactic 
center. By pushing these values to v c — 280 km/s and 
i?o = 9.4 kpc, we can barely accommodate a mass gap of 
8M = 2m e , but these choices are respectively 4ct and 3tr 
away from the mean values. 

One way to ease the tension would be to consider mod- 
els where the excited states are charged and thus able to 
decay into single electrons or positrons; this would lower 
the required mass gap to just 6M = m e . But this re- 
quires a great deal more model-building gymnastics than 
the case of a neutral excited state. A theoretically more 
appealing alternative is the case of three DM states where 
XoXo —5" XiX2 ; m which the mass gap SMqi between xo 
and xi is much smaller than SM 02 > 2m e 0, H3]. For 
example if SMqi = 0, then the effective 5M which de- 
termines the threshold velocity is half as large as in the 
models that produce two of the same excited state. This 
possibility can arise if the hidden sector gauge group is 
SU(2)xU(l) and gets completely broken, for example. 
However to see if this class of models can really have a 
larger rate would require solving the Schrodinger equa- 
tion in the three-state system, so while it is plausible that 
such asymmetric excitations could increase the rate, it is 
not proven by our analysis. 

Another way of getting around our no-go result is to 
assume the existence of a stable or metastable popula- 
tion of already excited states Xi) which need only be 



further excited to X2 through a small mass gap 8M\2, 
if 5Mq2 is assumed to be greater than 2m e so that the 
decay X2 — > Xo e + e~ is allowed. We refer to this as the 
inverted mass hierarchy scenario. This idea comes with 
new complications, since it is difficult to prevent the ex- 
citation process from occuring earlier in the history of 
the universe, to maintain the population of xi to the 
present day. However there is an intriguing possibility to 
overcome this by assuming the present DM particles are 
products of the late decay of a much heavier predeces- 
sor, and thus were initially relativistic, at a time when 
they would normally have been nonrelativistic had they 
been produced through conventional freeze-out [28] . The 
higher velocity suppresses the Sommerfeld enhancement 
at early times, when it would have the undesirable effect 
of depleting the metastable intermediate DM states. We 
hope to examine this scenario more carefully in future 
work. 

We have not tried to impose in detail the additional 
constraints on the model which arise if one would like it 
to also account for the PAMELA and Fermi/LAT excess 
electrons/positrons, although we did focus on TeV scale 
DM for that purpose. If the DM explanation of the high- 
energy leptons is not ruled out by upcoming analyses, 
this would be an interesting next step. It might also be 
worthwhile to investigate the effect of a nonspherically 
symmetric DM halo [48| on the rate of positron produc- 
tion. 

Acknowledgment: We thank Gil Holder for valuable 
discussions about the galactic parameters. 
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Appendix A: Born approximation cross section equation, we here give the quantum field theoretic ex- 

pression in the Born approximation. The spin-averaged 
For comparison with the nonrelativistic excitation squared matrix element is 
cross section from numerical solution of the Schrodinger 



l-^l 2 = 4 -9 4 \ n [ 2 ( s2 + " 2 ) + 4M +* - M + + M - ~ 6M 2 M 2 ] 



(t ,; 2 ) 2 

4 

(u-n 2 ) 2 
4 



[2(s 2 + t 2 ) + 4M|u - M% + Alt - 6Af 2 M 2 ] 
[4s 2 - 8M 2 s + 2M 2 s + 3(M 2 - Ml) 2 } 



(Al) 



where g is the couping constant, s, t, u are the Mandelstam variables, M± — Mi ± Mo, and Mi, Ma and fj, are the 
dark matter and gauge boson mass respectively. 
In the center-of-mass frame, 



a = E 2 crn , t = -2p 2 (l - cos(e)), u = -2p 2 (l + cos(0)) 

so the cross section is 

I f dVt cm 1 2 W.| A ^i2 



(A2) 



E cm v rel .1 47T 8lT E, 



9 



2nE^v re i 



s 2 - 16M 2 s + 4M 2 s + 6(M 2 - Ml) 2 n 2 + {p + p') 2 

v ..i i 7Z Zi\2 I 



cm-ret ^ pp' {[I 2 + p 2 + p' 2 ) [I 2 + (p - p') 2 

64^ 4 + 192^ 2 p 2 + 192^ 2 p' 2 + 64// 2 M 2 + 160p 4 + 160p' 4 + 192p 2 p' 2 + 32s 2 - 16M 4 + 16Mf - 96M 2 M 2 



16(2p 2 + 2p' 2 + M 2 +n 2 ) fi 2 + (p + p') 



{fi 2 + (p - p') 2 )(ji 2 + (p + p') 2 ) 
,) -)\ 



pp' n 2 + (p — p 

We rewrite the kinematic variables using the dimensionless quantity A, 



(A3) 



v rel = 2v t VK, p = MoVtVK, p = M v t VA - 1 



E r 



2M \/1 + v 2 A, M + = 2M Q [ 1 + If ] , M_ 



M 



(A4) 



Using »(<<1 and SM <C [i <C Mo, we obtain the approximate expression (fT3"|) for the cross section. 

I 



Appendix B: Escape velocity 

To be self-contained, we reproduce here the argument 
of ref. [lj| for the r-dependence of the escape velocity. It 
is important to notice that dark matter does not dom- 
inate the mass of the inner galaxy. To account for the 
gravitational effect of the baryons, one can use the fact 
that the circular velocity v c (the velocity of a test mass 
on a circular orbit around the galactic centre), is nearly 
constant out to radii or order r_2- We make the simpli- 
fying assumptions that there is no mass beyond r = r_2 
and that the density is spherically symmetric. 

The gravitational potential is 



-G 



M(r') 



dr'. 



where M(r), the mass within radius r, can be inferred 
from F = ma for a test particle of mass m: mv 2 /r = 
GmM{r)/r 2 . Hence M(r) = rv 2 /G for r < r_ 2 and 
M = r_ 2 v 2 /G for r > r_ 2 . The integral for $ becomes 



r-2 , 



(B2) 



The escape velocity is given by ^v 2 sc 



(Bl) eq. CE7J). 



$ = 0, leading to 
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